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ABSTRACT 



We report a possible change in the orbit parameters of the TrES-2 exoplanet. With a period of 2.470621 days, the TrES-2 exoplanet 
exhibits almost "grazing" transits 1 10.4 minutes duration as measured in 2006 by Holman and collaborators. We observed two transits 
of TrES-2 in 2008 using the 1.2m Oskar-Luhning telescope (OLT) of Hamburg observatory employing CCD photometry in an i-band 
and a near to R-band filter. A careful lightcurve analysis including a re-analysis of the 2006 observations shows that the current transit 
duration has shortened since 2006 by ^ 3.16 minutes. Although the new observations were taken in a different filter we argue that the 
observed change in transit duration time cannot be attributed to the treatment of limb darkening. If we assume the stellar and planetary 
radii to be constant, a change in orbit inclination is the most likely cause of this change in transit duration. 
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O ■ 1. Introduction 

c/^ , The study of transits is one of the most powerful methods in ex- 
■ ploring extrasolar planet properties. Although transits events are 
rare, viewing geometry dependent phenomena, they can provide 
H information about extrasolar planets and that is otherwise inac- 
^ ■ cessible. Space missions such as CoRoT or the Kepler mission 
' use the transit method in surveys of extrasolar planet systems. 
CO The transit light curve incorporates all relevant physical system 
parameters and, therefore, a highly accurate extrasolar planet 
transit light curve provides access to this information. The or- 
i/^ ■ bits of the planets in our solar-system as well as the Moon and 
" artificial Earth-orbiting satellites are known to undergo secular 
0^ changes in their orbit parameters. The reasons for these changes 
are the gravitational attraction of other ("third") bodies, devia- 
^ tions from spherical symmetry, air drag, non- gravitational forces 
• ^ and relativistic eflTects among others. Similarly, the orbits of ex- 
^ trasolar planets are expected to vary at some level, although, 
^ there is so far no evidence that the orbit or the physical parame- 
. . I ters of any exoplanet have changed. 

Changes in the orbit inclination / caused by a precession of 
the orbit plane are particular interest. Transiting planets are ideal 
for detecting these changes, especially when the transit is "graz- 
ing", i.e., when the planet eclipses only the polar regions of its 
host star. As a consequence, extrasolar planets with lower in- 
clination and larger impact values (but still producing transits) 
such as the cases of OGLE-56, TrES-2, TrES-3, and TrES-4 are 
particularly well suited to detecting these orbital changes. 

The TrES-2 exoplanet system was discovered in 2006 
by the TrES (Trans-atlantic Exoplanet Survey) project 
(iQ'Donovanet al . (2006)). Using t he Mt. Hopkins Observatory 
1.2m FWLO, iHolman et aP (l2007l) performed the first accurate 
analysis of TrES-2, using three TrES-2 light curves obtained 
in the fall of 2006. TrES-2 turned out to be a rather typical 
hot Jupiter exoplanet of mass Mp = 1.1 98M/ and radius 
Rp = 1.2227?/, orbiting a GOV star of mass = OMMsun 
and radius = l.003Rsun in a 2.47 day orbit. In this Letter, 



we present new transit observations of TrES-2 obtained on May 
and September of 2008, which suggest that the orbit inclination 
and hence the transit duration of TrES-2 have changed. 



2. Observations and data reduction 

We obser ved two transits of TrES -2 usin g the ephemeres sug- 
gested bv lQ'Donovanetal.1 (l2006h and bv lHolman elaP (l2007h 
from 

Tc(E) = 2,453,957.6358[if/Z)] -h E • (2.47063 days), us- 
ing the 1.2m Oskar Liihning telescope (OLT) at Hamburg 
Observatory. For the first observing run on 20 May 2008, we 
used SilKxlK CCD with a 5'x5' FOV, readout noise of 4.68 
e~, and gain of 1.12 e~ /ADU without filter (corresponding to 
a near R band). With a 30-second exposure time and 1x1 bin- 
ning, we achieved a time resolution of 1.13 minutes. For the sec- 
ond observing run on 18 September 2008, we used a new CCD 
camera with a 3Kx3K chip and lO'xlO' FOV, readout noise and 
gain were 16.37 e~ and 1.33 e~/ADU respectively, with an i- 
band filter. With this setup, we used a 60-second exposure time 
and a binning of 2x2, which provided an eflTective time resolu- 
tion of 1 . 17 minutes. During our observations, the airmass values 
ranged from 1.661 to 1.081 and from 1.0423 to 1.7176, during 
the first and second observing runs, respectively, and the seeing 
was typically 2.93" and 1.85", respectively, which are quite typ- 
ical of the Hamburg site. 

For the data reduction, we used Starlink and DaoPhot soft- 
ware routines, and the MATCH code. After applying bias sub- 
traction, dark correction, and flat fielding, we continued with 
our photometry. For TrES-2, we selected the aperture photom- 
etry mode using apertures centered on the target star, check 
stars, and sky background. Typical sky brightness values (the 
night sky is quite bright near Hamburg) for the first and sec- 
ond night was 325.6 and 270.2 ADUs, respectively, i.e., values 
at a level 1.5% and 1.25% of the star's flux, respectively. We 
used a total of 4 reference stars to test and calibrate the light 
curve (^71320+ 10215660-10214928-10220179-10219455). 
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Fig. 1. Observed TrE S-2 light curves, model fits and residuals (lower panel); upper three light curves were taken 1.2m FLWO 
(iHolman et alJ (l2QQ7b ). but fitted with our model code using Holman et al.'s (2007) limb darkening. The lower two light curves were 
taken with the 1.2m OLT at Hamburg Observatory, and fitted with the same model and the appropriate limb-darkening coefl&cients. 
The transit duration measured in 2008 is ~ 3 minutes shorter than that measured in 2006. 



Table 1. Relative Photometry data. 



HJD 



Relative Flux Uncertainty 



2454607.3512 
2454607.3528 
2454607.3536 
2454607.3543 
2454607.3551 
2454607.3559 



1.001980 
1.002020 
0.999300 
1.002150 
0.999920 
0.998690 



0.0015717 
0.0021024 
0.0014341 
0.0031083 
0.0030757 
0.0030854 



Note : Table 1 is available electronically. The time stamps refer to the 
Heliocentric Julian Date (HJD) at the end of each exposure. 



To esti mate the magnitude errors, we followed rHowell & EverettI 
(I200lh and used the expression 



CTmas = 1.8057 



(1) 



where p = ripi^il -\- -^)(Ns -\- No -\- -\- G^crj^), N, is the to- 
tal number of photons, Npix is the number of pixels, which was 
used to define A^*, hb is the total number of background pix- 
els, Ns is the number of total number of photons per pixel from 
background. No is the total number of dark electrons per pixel, 
A^^ is the total number of electrons per pixel due to the readout 
noise, G is the gain, and crj, is the Icr error of the A/D con- 
verter (~ 0.289). The factor of 1.805 7 converts the errors in flu x 
(electrons) into errors in magnitudes (iHowell & EverettI (l200lh ). 
we made no further changes to the adopted values of the mea- 
surement errors and continued with the light curve analysis with 
transit model fitting. Our final relative photometry is presented 
in Table [H which is available in its entirety in machine-readable 
form in the electronic version of this paper. 



3. Model analysis 

There are va rious ways of creating a physical tran sit model as 
described bv iMandel & Agoll (l2002h . IWinnI (l2008h . or lKippingI 
(l2Q08h . The main physical parameters of an exoplanet system 
are the host star mass and radius M^, R^, the planet mass and 
radius Mp and Rp, the orbit semimajor axis a orbits the orbit incli- 
nation i orbit, and the limb darkening law. There are also indirect 
or secondary parameters that can be computed once the main 
parameters are known, or which have no specific physical mean- 
ing or cannot otherwise be measured. These parameters include 
the central transit time Tc, which can be observed and used as a 
phase reference. The planetary period P must be determined by 
other means (i.e., not from the shape of the transit light curves), 
the eclipse duration and depths are functions of the orbit inclina- 
tion, the stellar and planetary sizes and the limb darkening law. 
In our modelling, we assume exactly spherical stars and planets 
with radii independent of the observing wavelength. We also as- 
sume circular orbits and - of course - the validity of Kepler's 
third law, which links the parameters a orbit, and Mp, since 
the period P is known. Usually one needs to fit all the main light 
curve parameters simultaneously, although, we clearly demand 
that the masses of star and planet, their radii, and the semi-major 
axis of the orb it remain constant for a ll ou r transit observations . 

Following Mandel & AgoU (l2002h and lHolman et al.l (l2007h . 
we assume a quadratic limb darkening law of 



J- = wi(l -Jd)- U2(l -fif, 



(2) 



where ji denotes the cosine of the angle between the surface nor- 
mal and the observer, and lo and Ip are the intensities at disk 
center and angle ji respectively; a complete list of the system 
parameters is given in Table O While in principle the limb dark- 
ening coefficients could be fitted by a transit light curve, limb 
darkening is indicated by the structure of the star's atmosphere 
and should not have to be adjusted to fit a specific transit light 
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Table 2. TrES-2 physical parameters (adopted from Holman et 
al.) 



TrES-2 system parameters 


Values 




0.98 Mo 


Mp 


1.198 My 


R, 


1.003 Ro 


Rp 


1.222 Rj 


Orbit Radius 


0.0367 At/ 


Period 


2.470621 days 



curve. Based on lClaretl (I2QQ4 and lSozzetti et alJ (l2QQ7h . we used 
the values ui = 0.318 and U2 = 0.295 for the i filter and ui = 
0.430 and U2 = 0.265 for the R filter. The relation between dura- 
tion and inclination is given by E q.[3] (.Charbonneau et al.. (.2006) : 
ISeager & Mallen-QrnelasI (l2003h ) 



D = — arcsin 

TT 




(3) 



where all parameters have the meanings explained above and D 
denotes the total transit duration. It is important to realize that 
the duration of a transit as defined by Eq.[3]is no t a function of 
the ass umed limb darkening coefficients (see also lKnutson et al.l 
(l2007h ). since the stellar and planetary radii are assumed to be 
wavelength independent; only the transit shape and depth, but 
not the transit duration depend on the limb dark coefficients. 

To determine the best fit model parameters, we used the 
statistic of 



X 



^ j,obs ■ 



j,mod 



(4) 



where Nf denotes the number of available data points, Fj^obs the 
flux, Fj^mod the model flux, and crj the error in Fj^obs all at time j. 
To ensure a consiste nt analysis, we re-analyzed the light curves 
already presented by lHolman et al.l (2001) using the appropriate 
limb darkening coefficients in addition to our own new transit 
observations. 

To determine the error in the derived fit parameters, we used 
the following bootstrap procedure, which has some similarities 
to a Markov chain error analysis. In step (1), we obtained the 
best-fit light curve minimizing Eq. |4l In step (2), we computed 
the best-fit model residuals FR defined by (FRj = F j^obs-F j^mod)- 
In step (3), we determined the second moment cr of the residual 
distribution, assuming a Gaussian distribution. In step (4), we 
generate a random light curve {RLj) using the model and the 
derived value of cr, and finally in step (5) determined a best-fit to 
this randomized light curve by repeating step (1). This procedure 
was repeated typically 1000 times and we recorded the variations 
in the fit parameters such as inclinations and durations. 

4. Results 

With the method described above we analyzed all available tran- 
sit observations (cf., Eig. [T]); the results of our fits are summa- 
rized in Table [3l where we also indicate the darkening coeffi- 
cients used in the analysis (upper five entries). Eor the Tabled 
we used three sets of limb dark coefficients. Set-1 for Holman's 
data (Sl:wi = 0.22, U2 = 0.32), Set-2 for OLT's data in R filter 
(S2\ui = 0.430, U2 = 0.265), and finally, Set-3 for OLT's data in 



Table 3. Duration, inclination,;^^ values and limb darkening co- 
efficients from five light curve fits; units for duration and errors 
are minutes and for inclination and errors are degrees. 



T/- tim^ JTy\ 

ic 1 11 lit ynjLJ\ 


LJUl UllUii 


Eji i Ui 3 


InLLLn. 


Eji 1 (Jl i3 


^ vUiUt: 


T JIT 


3989.7529 + 0.0069 


110.308 


0.432 


83.59 


0.019 


432.1 


SI 


3994.6939 + 0.0066 


109.230 


0.448 


83.56 


0.019 


296.8 


SI 


4041.6358 + 0.0070 


109.025 


0.430 


83.55 


0.019 


290.6 


s\ 


4607.4036 + 0.0072 


106.620 


0.883 


83.44 


0.036 


179.1 


S2 


4728.4740 + 0.0071 


106.112 


0.870 


83.43 


0.036 


190.1 


53 


3989.75286 


113.275 




83.74 




443.0 


S2 


3989.75286 


111.632 




83.66 




435.6 


53 


3994.69393 


111.450 




83.65 




305.5 


S2 


3994.69393 


110.195 




83.60 




300.1 


53 


4041.63579 


111.790 




83.67 




296.6 


52 


4041.63579 


110.197 




83.60 




292.4 


53 


4607.40356 


104.317 




83.35 




180.0 


51 


4607.40356 


105.060 




83.38 




179.1 


53 


4728.47400 


104.890 




83.37 




195.7 


51 


4728.47400 


107.100 




83.47 




192.6 


52 



i filter (S3:wi = 0.318, U2 = 0.295). In the first part of table, we 
listed the best-fit results using the most suitable set of limb dark 
coefficients. In the second part, we used all the possible com- 
binations of limb dark coefficients to show that duration is not 
strongly aff'ected by limb darkening. Holman's light curves have 
433, 301, and 299 data points, and ours have 169 and 193 data 
points, respectively. With our new observations, we now have 
precise measurements of Tc for five transits, so we are able to 
refine the ephemeres equation. A least squares fit to these five 
transit events results in the expression 
Tc{E) = 2, 453, 957.63403 [HJD] -h E • (2.4706265 days). 




Fig. 2. Eit improvement Ax^ for 1000 Monte Carlo realisations 
when going from models of fixed inclination to models of free 
inclination (see text for details). 

In our basic analysis, we assumed the period, masses, radii, 
and semi-major axis of the orbit to have the values derived by 
[Holman et al. ( 2007). With this approach, we found that the OLT 
data taken in 2008 appear to prefer diff'erent inclinations than the 
ELWO data taken in 2006 (cf.. Table O. To assess the signifi- 
cance of these inclination diff'erences, we carried out the follow- 
ing analysis. We assumed that the five transit light curves (cf., 
Eig. [T]) can be described by variable stellar and planetary radii, 
variable inclination, and variable central transit times (model A; 
20 fit parameters), and generated 1000 realizations of this light 
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curve with the method described above. In model B (10 fit pa- 
rameters), we fixed the stellar and planetary radii to a common 
value, but still allowed the inclinations and central transit times 
to vary for each data set, while in model C with (5 fit parame- 
ters), we also kept all inclinations constant at the best-fit model 
value and allowed only the transit times Tc to vary. The gener- 
ated random light curves were then fitted (via minimisation) 
by models (A), (B), and (C), and the errors were computed as 
described above; clearly, the closest fits are expected for model 
(A), and the poorest fits for model (C), which is, however, the 
"canonical" model. The essential question is whether the fit im- 
provement between models C, B, and A is significant or not ? 

We first consider the fit improvement between models (C) 
and (B); for each generated light curve, we can compute the 
appropriate A;^^, and the distribution of these /^x^ residuals is 
shown in Fig. [21 A typical value of /^x^ is clearly ^ 200, and 
an F-test indicates that the fit improvement between model (C) 
and model (B) is highly significant; carrying out the same exer- 
cise for going from model (B) to model (A), we find any addi- 
tional fit improvement is insignificant. This does of course meet 
our physical expectations, since the sizes of the star and planet 
are not expected to change, while the orbit inclinations might 
change. 

The bootstapped error distribution of the derived inclinations 
in model (B) at the five epochs is shown in Fig.[3l which shows 
the distributions to have very little overlap. We thus conclude 
that from a statistical point of view the change in inclination 
observed in the 2008 data is very significant. 



0.25 




83.3 83,4 83.5 83,6 83.7 

Intlination ["] 

Fig. 3. Inclination error distribution derived from 1000 Monte 
Carlo bootstraps: OLT light curves (left) and FLWO light curves 
(right). 



Can the apparent change in inclination be caused by errors 
in the analysis? We investigated possible non-physical causes 
of the derived inclination change. Our own OLT observations 
we re not carried out w ith the same filter as the observations 
by iHolman et al.l (l2007l) . Since our fit approach explicitly as- 
sumes that the stellar radius is independent of wavelength, the 
eclipse duration defined by Eq.[3]depends only on both the sizes 
of the host star and planet and the orbit inclination, but not on 
the shape of the transit light curve, which depends on the cho- 
sen limb darkening law. Since we keep the stellar and planetary 
radius fixed in our fitting scheme (model B), a change in incli- 
nation is the only way to produce a change in eclipse duration. 
However, since the ontire observed light curve is fitted to an as- 
sumed model light curve, we can obtain an incorrect fit and thus 



an incorrect inclination and duration if we assume an inappro- 
priate limb darkening law. 

Since the Hamburg OLT observations were carried out in 
"bluer" filters (than the filter which Holman used) with larger 
contrasts between disk center and limb, we assessed the influ- 
ence of limb darkening by analyzing all data for all limb darken- 
ing coefficients, assuming a quadratic limb darkening law spec- 
ified by Eq. [2l the results of these fits are also given in Table 
[3] (lower ten entries). As is clear from Table O the derived du- 
rations depend on the assumed limb darkening coefficients, al- 
though, the inclination discrepancy between the 2006 and 2008 
data cannot be attributed to limb darkening. We further note that 
the nominal limb darkening coeflftcients provide the best fits in 
all cases. 

Finally, we asked ourselves which sets of limb darkening 
model coefficients would provide the best-fit light curve to our 
new OLT data, keeping the transit duration fixed at the val- 
ues observed in 2006. With this fit setup we found poorer fits 
ix^ = 204.58) with unrealistic limb darkening coefficients of 
ui = 0.4258 an d U2 = 0.2254 , corresponding to Te/f ~ 4800^ 
for epoch-312 (IClaretl (l2004l) ). too low given the spectral type of 
TrES-2's host star. We repeated the experiment for all combina- 
tions of limb darkening coefficients in the range 0.00 < -h ^2 < 
1.00, and found a best-fit model with;^^ = 182.0, although, the 
corresponding coefl&cients of = 1.19 and U2 = -0.80 do not 
relate to any significant physical meaning. We therefore con- 
clude that the discrepancy between the derived durations of the 
2006 and 2008 transits cannot be explained by the assumed limb 
darkening laws. 

Another possible source of error is an incorrect determina- 
tion of the out-of-transit flux. All assumed model light curves as- 
sume an out-of-transit flux of exactly unity. The data points are 
normalized to unity by the available out-of-transit data points, 
and therefore the accuracy of this normalization depends on the 
number of such out-of-transit data points and their respective er- 
rors; we note that in the extreme case of zero out-of-transit data 
points no normalization to unity can be carried out and therefore 
no meaningful analysis. If the data normalization were for some 
reason too low, the observed transit depth would be deeper and 
hence the derived inclination higher (and duration longer); what 
must then concern us is a data normalization that is too high, re- 
ducing the eclipse depth and hence reducing inclination and du- 
ration. To study the influence of errors in the data normalization, 
we introduced a multiplicative amplitude factor to the model, 
which should nominally be unity but which was allowed to vary 
to obtain a best fit. In all cases, we found that these amplitudes 
were close to unity (at the level of a few times 10""^) without 
any substantial change in transit duration. We thus conclude that 
normalization errors do not significantly aff'ect our results. 

Finally, following * Holman et al.l (12007), our modelling ap- 
proach assumes a spotless star. If the true star were spotted, 
our modelling approach and thus the inferred model parameters 
would be incorrect. Clearly, if one were to place spots exactly at 
those limb positions where the planets enters and leaves the stel- 
lar disk, the apparent transit duration would change; at the same 
time the limb darkening and transit depth would change and it 
is unclear, what the resulting fit parameters of such an inappro- 
priate model would be. Since we, first, have two transit obser- 
vations yielding consistent fit results and, sec ond, the vsinji) of 
TrES-2's host star is quite low (< 2.0 /:m/>s'gc. [b' Donovan et al.l 
(2006)) proposed that it was an inactive star and we conclude 
that star spots are an unlikely explanation of the observed transit 
duration shortening. 
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Fig. 4. Epoch versus inclination: a linear fit is clearly sufficient 
to describe the available data; if the inclination continues to de- 
crease, the minimum inclination threshold (imin = 81.52^) will 
be reached around epoch -4179. 
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If we accept that the formally measured inclination changes 
are real, we can study (in Fig.© the variation in inclination ver- 
sus time derived from the presently available data; clearly, a lin- 
ear fit is sufficient for a description. In every transit systems, 
two theoretical thresholds of minimum inclination angles exist, 
a first threshold, where the planet is never fully in front of the 
stellar disk, and a second threshold, where no transit occurs at 
all. For the TrES-2 system parameters listed in Table [2l the first 
inclination threshold occurs for imin,i = 83.417^, and the second 
threshold for imm,2 = 81.52^. From our linear fit shown in Fig. HI 
the first threshold is reached around Epoch - 325 corresponding 
to October 2008 with an eclipse duration ofD= 105.92 minutes. 



5. Conclusions 

Because its orbital properties of TrES-2 is a particularly interest- 
ing exoplanet. The transits of TrES-2 are almost "grazing" and 
are therefore quite sensitive to inclination changes. By compar- 
ing our new 2008 observations of TrES-2 with data obtained in 
2006, we have detected an inclination change of 0.1 degress in a 
little under two years, corresponding to a transit duration change 
of about ~ 3.16 minutes. We have investigated the influences 
of limb darkening and the light curve normalization on the de- 
rived inclinations and found that the observed duration changes 
cannot be attributed to incorrect light curve modelling. We had 
therefor to search for a physical cause of the observed changes. 
Since stellar and planetary radii are unlikely to have changed 
and an explanation related to star spots appears unlikely given 
the activity status of TrES-2's host star, a change in the orbital 
inclination is the only remaining explanation. A third body in 
the form of an outer planet, which possibly does not cause any 
transit eclipses, would provide the most straightforward expla- 
nation. We emphasize that it appears extremely worthwhile to 
continue transit monitoring of TrES-2 in the 2008 - 2010 time 
frame to study its inclination changes as well as to search for fur- 
ther companions of TrES-2 in radial velocity data. In the Kepler 
field, high precision photometry of hundreds of transits of TrES- 
2 will be obtained, which should allow a precise determination 
of all orbit changes as well as a refinement of the stellar and 
planetary parameters of the TrES-2 exoplanet system. 
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